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We consider correlated Levy walks on a class of two- and three-dimensional deterministic self- 
similar structures, with correlation between steps induced by the geometrical distribution of regions, 
featuring different diffusion properties. We introduce a geometric parameter a, playing a role 
analogous to the exponent characterizing the step-length distribution in random systems. By a 
single-long jump approximation, we analytically determine the long-time asymptotic behaviour of 
the moments of the probability distribution, as a function of a and of the dynamic exponent z 
associated to the scaling length of the process. We show that our scaling analysis also applies to 
experimentally relevant quantities such as escape-time and transmission probabilities. Extensive 
numerical simulations corroborate our results which, in general, are different from those pertaining 
to uncorrelated Levy-walks models. 



I. INTRODUCTION 

Levy-like motions represent an important family of 
random motions, generated by a stochastic process with 
stationary and independent increments. Brownian mo- 
tion constitutes a particular case of this family, but its 
mathematical tractability and its remarkable statistical 
properties have led it to become the principal model of 
random motion. However, different types of Levy mo- 
tions have been extensively investigated, both theoret- 
ically and experimentally, as they have been found to 
be ubiquitous in nature: in biology [11-3 j chaotic dy- 
namics d, II] , economics , search strategies [1, , the 
Levy motion that have been observed are characterized 
by increments of arbitrary length I, with a step length 
distribution featuring an algebric tail ^ Such a 

distribution is said to be heavy-tailed and has a diverging 
variance for a < 2. 

An interesting experimental situation where Levy mo- 
tion can be detected is diffusion in heterogeneous and 
porous materials, composed of two o more types of re- 
gions with different diffusion properties. In that case, 
the motion of particles consists of a sequence of scatter- 
ing events occurring in the hard-scattering part of the 
material, followed by long jumps performed at almost 
constant velocity in the non-scattering regions. If the ma- 
terial is very heterogeneous on all scales, the step length 
results to be heavy-tail Levy-distributed. This descrip- 
tion applies to transmission of light through clouds [l^l , 
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tracer transport in heterogeneous aquifers 
ular diffusion at low pressure in porous media 
is dominated by collision with pore walls, with ballistic 
motion inside the large pores [12| -, as well as in frac- 
tured and heterogeneous porous media [l^ . In addition, 
recent experiments on new disordered optical materials, 
the Levy Glasses, paved the way to the engineering of 
Levy-distributed step-lengths [ij. This phenomenology 
is often modeled in terms of annealed Levy walks, where 
step- lengths are uncorrelated [H, [11] . 

However, a key signature of these processes is that the 



steps are in principle not independent, as they are corre- 
lated by their mutual positions in the sample. A walker 
that has just traversed a large hole has a high proba- 
bility of being backscattered at the following step and 
thus to perform a jump of roughly the same length. The 
step length distribution represents therefore a " quenched 
disorder" , a standard definition in statistical mechan- 
ics of disordered systems [l^]- Now, while the case of 
Levy walks with uncorrelated jumps is well understood 
Tsl . i^[iM[i^ I the correlation effects, which are expected 



to exert a deep influence on the diffusion properties, are 
still to be characterised, and the upper critical dimension 
above which a description in terms of uncorrelated Levy 
walks proves sufficient is still under debate [20l-|23j. 

A first step in this direction has been taken, and to 
this end, quenched Levy processes have been studied on 
one-dimensional systems [23,[2^. There, the effects of 
geometry-induced step-length correlation on the asymp- 
totic behavior of the mean square displacement has been 
investigated in detail. More recently, again for one- 
dimensional systems, different aspects regarding the scal- 
ing properties of random- walk distributions, the relations 
between the dynamical exponents and the possible aver- 
age procedures have been discussed in a common frame- 
work. This evidenced the strong effects of correlations 
also on quantities averaged over all the starting sites (26l - 
[28j . Obviously, one-dimensional models represent simpli- 
fied systems and may not compare quantitatively to real 
experiments, so far conducted on three dimensional sys- 
tems. On the other hand, they have proven amenable to 
an exact analytic solution for the dynamics. It is now 
important to keep track of what properties observed in 
one dimensional systems can be extended to higher di- 
mensions. The key point is to understand the effects 
of step-length correlations and of averages over starting 
points, and what is needed is a proper model for the 
heterogeneous geometrical pattern where the correlated 
Levy motion takes place. 

In this paper we model the transport process as a ran- 
dom walk on scale-invariant structures — namely gener- 
alized two- and and three-dimensional Sierpinski carpets 
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— whose solid and empty regions represent portions of 
material characterized by different transmission proper- 
ties. In order to reproduce a realistic propagation, our 
process consists of a simple diffusive random walk on the 
solid regions of the structure, and of a ballistic motion 
across its empty regions. Therefore, while the distribu- 
tion of allowed jump lengths is scale-invariant, any point 
in the two- or three-dimensional region of space occupied 
by the structure is accessible to the walker. This is at 
variance with the well-known and studied case of stan- 
dard diffusion on self-similar structures, where the empty 
regions are not accessible [l^ , and the motion takes place 
on a truly fractal space. From a more rigorous point of 
view the ballistic propagation of the walker across the 
"empty" regions of fractal structure can be equivalently 
obtained by imposing a scale invariant pattern of persis- 
tence on the sites of a regular Euclidean lattice [3^ ■ 

It should be remarked that the available experimen- 
tal realizations are characterized by a strong amount 
of disorder (lol - [l3 |. which might call into question our 
modeling of the pattern of inhomogencity in terms of de- 
terministic, self-similar structures. However, the regular 
distribution of the inhomogeneities allows us to take into 
account several crucial points in a more controlled way. 
We are able to give an estimate of the effects of the tail 
in the step length distribution in terms of the geomet- 
rical parameters characterizing the Levy quasicrystals, 
and we can tunc it arbitrarily to study its effects; we 
can correctly control the effects of the averages over the 
starting sites, which is known to have a deep influence the 
asymptotic properties of the mean square displacement 
in one-dimensional structures; most importantly, we are 
able to take a first step towards the comprehension of 
correlated Levy walks in higher dimensions. 

One of our main results is that, for this type of "topo- 
logical" correlations, quenched Levy walks lead to a set 
of exponents for the asymptotic behaviors of the mean 
square displacements and transmission which differ from 
the uncorrelatcd case [20l42^ . 

In detail, in this paper, we tackle the problem by verify- 
ing the scaling hypothesis for the probability distribution 
Pj (r, t) , namely the probability for the Levy walker to be 
at a distance r from its starting site j at time t [IH and 
by evaluating the dynamical exponent z associated to the 
growth of the scaling length £{t) ^ i^/^ of the process. 
This is a function of a tunable geometrical parameter, 
a, describing the self similar pattern of empty regions in 
the Levy quasicrystals, and playing a role analogous to 
the exponent characterizing the step-length distribution 
in random systems. Then, by making the reasonable and 
well verified "single long jump" hypothesis (27| . we esti- 
mate the tails of P(r, <), i.e. the average of Pj (r, overall 
the possible starting sites. We derive an analytic expres- 
sion for the asymptotic behaviors of the moments (r^it)) 
of P{r, t), evidencing the presence of strongly anomalous 
diffusion [13, i.e. {rP{t)) ^ Ce{t)P. In the attempt of 
making a direct contact with experiments, we also de- 
rive the scaling properties of exit times, as well as of 
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FIG. 1: Third generation of a SC with tiu = 3 and rir = 8. 
Light squares represents actual sites of the structure. Dark 
squares belong to empty regions. The "exploded view" in the 
lower left corner highlights the fact that this structure is made 
of rir — 8 SCs of generation g — 2, and so on. 



the time-resolved transmission probability and transmis- 
sion profiles through a slab of thickness L. Extensive 
numerical simulations are in very good agreement with 
the predicted behaviors. 

The paper is organized as follows: in the following sec- 
tion we describe our Levy quasicrystals based on general- 
ized Sierpinski carpets (SC). In section Hill we introduce 
the relevant physical quantities that we will study in this 
paper. In Section IIVI we discuss the scaling hypothesis 
and the single-long-jump approximation which allows us 
to evaluate the momenta of the distribution P{r, t) evi- 
dencing the presence of strongly anomalous diffusion ■ 
In the last part of the paper we present extensive nu- 
merical simulations proving the reliability of our scaling 
hypothesis and single long jump approach. In particular, 
we evaluate the dynamical exponent z as a function of 
the dimensionality, of the dynamics and of a, a simple 
parameter describing the topology of the structure. For a 
slab of thickness L we evaluate exit times, time resolved 
transmissions and transmission profiles, such a results 
could be useful for a comparison with experiment 
Section IVll contains our conclusions and perspectives. 



II. STRUCTURES 

In order to to realize a Levy walk on a quenched struc- 
ture we consider a regular lattice in which empty regions 
have been created by removing a subset of sites. We then 
consider a random walker that can jump across empty re- 
gions. This is actually a Levy process provided that the 
linear size of such regions — and hence the length of the 
ballistic jumps — is Levy-distributed. Specifically, it is a 
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Levy walk if all jumps arc performed at constant veloc- 
ity, and hence completed in a time proportional to their 
length. The Cantor lattices considered in Ref. [l^l, are 
deterministic one-dimensional structures exhibiting the 
desired properties. The natural two-dimensional gener- 
alization of such structures is provided by the so-called 
Sierpinski carpets discussed in Ref. [s^ . 

The self-similar character of a SC can be illustrated by 
observing that several identical objects can be assembled 
to give a larger SC. One way to describe this constructive 
procedure is the following: arrange n.^ SC of generation 
g in a square array, and subsequently remove some, so 
that only rir < of them remain. The resulting object 
is a SC of generation 5 + 1. If the same procedure is re- 
peated using SCs of generation g a.s building blocks, 
a structure of generation g -\- 2 is obtained. Figure [T] il- 
lustrates this algorithm in the case of the standard two 
dimensional SC, where = 3, = 8 and the struc- 
ture of generation g = is a simple lattice site. At each 
generation a 3 x 3 array of carpets is laid out, and the 
central structure is removed, leaving a square empty re- 
gion. Note that the linear size of the largest empty region 
in the structure grows by a factor Till the generation is 
increased. A different structure, corresponding to = 4 
and rir = 12, can be obtained by laying out a 4 x 4 ar- 
ray of building blocks, and subsequently removing the 
four innermost structures. In general a SC of genera- 
tion g features n^r^^^ empty regions of linear n,^, with 
fc = 0,l,-- - ,(7—1. Ifgis sufficiently large, the linear size 
of the empty regions spans several orders of magnitude. 
For instance the 3rd generation structure in Fig [T] fea- 
tures one empty region of linear size 9, 8 empty regions 
of linear size 3 and 64 empty regions of linear size 1 . 

One of the parameters characterizing a SC is its frac- 
tal dimension, i.e. the average number of sites per unit 
surface, which is has a simple expression in terms of the 
constructive parameters 



a Levy walk in d dimensions, whose p-th moment 



log rir 
log nu 



(1) 



A more relevant characterization for our purposes can 
be given in terms of the probability for the length of the 
jumps that can be taken at a random site of the structure. 
For each of the possible jump lengths n,^ this probability 
is given by density of sites on the boundary of empty 
regions of linear size n^. It is easy to prove that the 
moments of such probability 
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converge only if p < a, where 

log tlr 



a 
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1. 



(2) 



(3) 



This parameter has exactly the same role as the exponent 
a characterizing the distribution of the step lengths s of 
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(4) 



converges only if p < a. 

For the Sierpinski carpets discussed in Ref. [s^ the 
threshold in Eq. ^ is connected to the fractal dimension 
of the structure, df = a + 1 < 2. Since a < 1 the size 
of the average empty region, i.e. the first moment of the 
distribution, is infinite. 

The two dimensional counterparts of the Cantor- 
Smith- Volterra sets discussed e.g. in Ref. [l^l are hyerar- 
chical structures similar to Sierpinski carpets for which 
a > 1 and df = 2. We refer to these structures as fat 
carpets (FC). 

As we discuss above, a SC of generation g can be seen 
as a ri„ X ri„ square tiling whose elements have the same 
size as SC of generation g — 1. Some of these elements 
are actual SCs, while others are left empty. The tiling in 
a FC is sligthly more complex, in that it requires three 
kinds of tiles, as it is illustrated in Fig. O This comes 
about because, in order to attain a > 1, the largest empty 
tiles at a given generation must grow less than the entire 
structure as the generation is increased. More specifi- 
cally, a FC of generation g + 1 contains i^^ FCs of 
generation g, whose side is \ {vr ~ 1)^ empty square 



tiles of side l'^^ and 2vr{vr — 1) rectangular L^^' by 



(s) 



-(9) 



FIG. 2: Third generation of a FC with rir = 9 and n„ = 2. 
The "exploded view" in the lower left corner illustrates the 
tiling pattern underlying the structure. Three different types 
of tiles are necessary to build the structure: Ur = Vr = ^ FCs 
of generation g = 2, whose side comprises 19 sites, [ur — 1)^ = 
4 empty square tiles of side — A and 2{ur — l)vr = 12 
rectangular padding tiles, highlighted in a lighter shading. 
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tiles. If L^Q^ ~ and L^'' = n^, with n„ < i/^, we get 



v — h 



-h 



(5) 



so that the fraction of sites contained in the rectangular 
padding tiles decreases as nf^v~^ . Thus, for sufficiently 
large generations, the contribution of the padding tiles 
becomes negligible and the total number of sites in the 
structure is proportional n-j. . This means that the thresh- 
old for finite moments in the probability for the length of 
the jump at a random site of the FC is once again given 
by the quantity in Eq. ([S])- Since Vr ~ ^fn^ > n„, a > 1 
and the average size of the empty regions is finite. All 
of the above discussion can be easily generalized to arbi- 
trary dimension d. In fact, in Section ITVl the single-long- 
jump picture is illustrated for a generic dimension d. As 
to the simulations, we mainly focus on two-dimensional 
structures, but we also present some results for three- 
dimensional sponge-like structures. 

Let us now get back to the diffusion dynamics of the 
random walk. This is based on the standard dynamics 
on a regular square lattice: one direction out of the four 
possible ones pointing to the adjacent sites is randomly 
chosen. If a site is found in that direction at a distance of 
a single lattice spacing, then the jump is taken, and the 
time is increased by one unit. If the same jump would 
land inside an empty region, then the random walk tra- 
verses the entire region and lands at a site belonging to a 
different boundary from that containing the starting site. 
Fig. [3] illustrates the two possible schemes for the choice 
of the arrival site we consider in this work. In the simpler 
scheme (left panel) the walker lands at the site facing 
its current position {head-on dynamics). In the second 
scheme (right panel) the landing site is randomly chosen 
among those belonging to any of the sides of the empty 
region except the one currently hosting the random walk 
(fan-out dynamics). In any case, since the random walk 
moves at constant velocity, the time required for a jump 
equals its length. 

Note that the jump probability of these two dynamics 
are different at any distance scale. Therefore we expect 
that parameters usually exhibiting universal character 
might be different for the two cases. In fact, our sim- 
ulations evidence that the dynamical exponent govern- 
ing the growth of the characteristic length of the process 
seems to depend on the dynamics, on the dimensionality 
of the system and on the geometric parameter a. Con- 
versely, a very robust feature we observe is that anoma- 
lous diffusion takes place strictly for a < 1 for every case 
considered here. This is at variance with the case of un- 
corrclated Levy walks, where a < 1 and 1 < a < 2 corre- 
spond to purely ballistic propagation and superdiffusion, 
respectively. 
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FIG. 3: Possible jumps of a random walk initially at a site on 
the boundary of an empty region (black circle). The random 
walk can attempt a jump in one of the four possible directions, 
with equal probabilities. If the chosen direction points to one 
of the three neighbouring sites (black triangles), the jump is 
taken and the time is increased by one unit. If the chosen 
direction points to the empty region, the random traverses 
it. The two panels illustrate the two possible schemes for 
the choice of the arrival site (black pentagram). Left: the 
arrival site is the one facing the current position of the walker 
{head-on dynamics). Right: the arrival site is equiprobably 
chosen in the boundaries of the empty region not containing 
the current position of the random walk {fan-out dynamics). 
In any case, the time is increased by the length of the jump, 
as measured in units of the lattice spacing. 



III. PHYSICAL QUANTITIES 

The diffusion process on the structures under concern 
can be characterized by several quantities. One of the 
most natural is the probability Pj (r, t) that at time t the 
random walk is at a distance r from the lattice site j it 
started from. The average of this probability over all the 
lattice sites, 



(6) 



where N is the number of lattice sites, is a meaningful 
quantity as well, in view of the spatial inhomogeneity 
of the SCs and FCs H^]. As we will discuss shortly, it 
highlights the anomalies in the behaviour of the diffusive 
process. These probability distributions are completely 
defined by their moments, 

{rP{t))= [ drrP Pj{r,t), {rP{t))^ j drr^ P{r,t) (7) 



A further quantity worth analyzing is the escape time 
probability, namely the probability £{L,t) that a walker 
starting at a random site of a structure of linear size L 
escapes from it at time t. 

An even more direct contact with the experiment is 
provided by the conditional first-passage time Tjk{L,t), 
i.e. the probability that a particle injected at a site j 
belonging to one boundary of the structure hits a site k 
belonging to the opposite boundary at time t before go- 
ing back to the boundary containing j, which is therefore 
absorbing. This means that there is a (large) probabil- 
ity that the injected particles are backscattered by the 
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structure. As to the direction parallel to the entrance 
and exit boundaries, we consider two possible schemes. 
In the first, the random walk is simply reflected by the 
boundaries of the structure. In the second the environ- 
ment seen by the random walk is periodic, i.e. it consists 
of an infinite stack of identical structures of some fixed 
generation. 

The conditional first-passage time gives access to ex- 
perimentally relevant quantities. For instance, it is 
clearly related, through integration over time, to the 
transmission profiles shown in Ref. (l^ . i.e. the average 
number of processes transmitted through the sample at 
a transverse distance r from the injection point. Specifi- 
cally, on a 2D system of linear size L, these are obtained 
as 

1 f°° 

%i'^^L)^-Y. I dtT,,+riL,t) (8) 

Likewise, averaging over entrance sites and summing 
over exit sites gives the the time-resolved transmission 
probability, 

r(L,t) = i^T,fe(i,i). (9) 

Finally, integrating this quantity over time gives the to- 
tal transmitted intensity, i.e. the fraction of processes 
transmitted through the sample 

/>oo 

I{L)= / dtT{L,t). (10) 
Jo 

IV. SCALING ANALYSIS AND SINGLE LONG 
JUMP ANSATZ 

Several of the quantities discussed in the previous Sec- 
tion depend on time and on a spatial variable. The latter 
can be the walk distance r in the case of the local and 
average probabilities Pj{r,t) and P{r,t), or the system 
size L in the case of the escape probability. A generic 
function of space and time is expected to exhibit the fol- 
lowing scaling structure 

where the asymptotic behaviour of the characteristic 
length is 

e{t) - (12) 

Notice that the prefactor in r.h.s. of Eq. (fTTj) depends 
only on time. An alternative, equivalent form of this 
equation is also possible, featuring a prefactor depending 
on the spatial variable alone, r"'^ , and a scaling function 
simply related to /. 

The exponent 7 in the prefactor can be determined 
from known features of the quantity under concern. For 



instance, the expected dependence of the return prob- 
ability on the spectral dimension dg of the structure, 
P(0,i) - t-'^"/^ [H, [H, entails that 7 = 4/2 for 
the local probability Pj{r,t). Furthermore, the time- 
independence of the norm of Pj{r,t) yields z = 2d/ds, 
where d is the dimensionality of the system. 

As to the function f{x), one expects it to decay very 
rapidly unless, at time t, the random walk was able to 
take jumps much longer than £{t). The rapid decay is 
typical of local quantities, whereas average quantities 
usually exhibit an anomalous behaviour. We illustrate 
the origin of the slow decay by focusing on Pj (r, t) and 
its average P{r,t), defined in Eq. (|6]). A slow decay of 
f{x) is expected when the random walk is able to take 
arbitrarily large jumps at arbitrarily small times. This 
is precisely the case of P{r,t), since the sum in Eq. © 
includes starting sites in the vicinity of empty regions of 
arbitrary size. 

The situation is somewhat different for the local prob- 
ability Pj (r, t) . It is indeed true that the starting site can 
be close to an empty region of very large linear size, say 
L, so that a jump of length L can occur at relatively early 
times. However, this merely affects the transient part of 
the dynamics. After a suitably long time £{t) ^ L, the 
regular behaviour is recovered, since empty regions much 
wider than £{t) are far removed from the position of the 
random walk. For a similar reason we expect the long- 
time asymptotic behavior of Pj (r, t) to be independent of 
j. Indeed for £{t) much larger than the distance between 
sites i and j we expect that Pj{r,t) ~ Pi{r,t). 

The fast or slow decay of f{x) is mirrored by the 
asymptotic behaviour of the moments defined in Eq. (O . 
The rapidly decaying scaling function in the local prob- 
ability indeed yields 

(rAt))- i/{r^it))-m (13) 

independent of the starting point. 

The expected more complex behaviour of the average 
moments (r''(t)) when slow decay is pres ent can be stud- 
ied in the single-long-jump picture j26l428| . This corre- 
sponds to assuming that the average probability can be 
split into two contributions, a regular, rapidly decaying 
part P'-''-' (r, t) and an anomalous tail P^'^'f (r, t). The latter 
is assumed to be generated by a single jump, much larger 
than the characteristic length, and it is expected to give 
rise to a significant contribution in the region r ^ i{t) 
where the former is negligible. 

In order to estimate the anomalous contribution, we 
have to determine the probability that the random walk 
takes a jump of size £ ^ £{t) at time t. This is pro- 
portional to the density of regions of linear size £ times 
the number of sites from which a random walk can reach 
one of such regions in a time t. Recalling that the aver- 
age distance covered by the walker in a time t is given 
by £{t), the sites accessible to a single region of linear 
size £ are enclosed in a region whose volume is propor- 
tional to £{t) £'^~^ . The density of sites corresponding to 
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a volume of linear size £{t) is £{t)'^'^'^, where d( is the 
fractal dimension discussed in Section |IT1 Observing that 
the empty regions in the self similar structures thereby 



discussed have discrete linear sizes i 



and discrete 



densities 



we estimate the anomalous contribution 



to the average probability as 



5[r-mm{ntt)] (14) 



where the delta function takes care of empty regions 
whose linear size is larger than the longest possible dis- 
tance (ballistically) covered by the random walker in a 
time t, and 



1 



a if a < 1 
1 if a > 1 



(15) 



We can now analyze the generic moment of the average 
probability. According to its definition, the regular part 
of such probability results in a regular contribution, like 
in Eq. i^: (rP)(') - £{t)P . The contribution from P'-''^ 
can be evaluated as 



rP\(a) 



^(*)<"fi<t 



(16) 



Performing the summations and joining the regular and 
anomalous contributions, we obtain 



rP) ^Cit- + C2 t- 



(17) 



where the Cj's are time- independent coefficients and we 
made use of Eq. (|12p . Equation P7)) shows that the 
asymptotic behaviour of the generic moment of the aver- 
age probability is actually more complex than that given 
in Eq. 1131 i. e. the system exhibits strongly anomalous 
diffusion [32|- More specifically on thin carpets (a < 1) 



lip < a 
ifp > q; 



whereas on fat carpets (a > 1) 



it 



if p < 2a 
if p > 2a 



(18) 



(19) 



Let us now turn to the scaling properties of the es- 
cape time probability. The scaling hypothesis states that 
£{L, t) can depend on the system size L only through the 
ratio L/£{t) i.e. 

£{L,t) = L\f{L/£{t)). (20) 
Imposing the normalization 

1= / £{L,t)dt= [ L^f{L/l{t))dt, (21) 



we get r] = —z, where we once again made use of Eq. (|12l) . 

Similar considerations apply to the transmission- 
related quantities introduced in Section IIIII The time- 
resolved transmission probability, Eq. ^ , is expected to 
obey a scaling law of the same form as Eq. (pO)) . 



Tit,L) = L^f{L/em 



(22) 



although with different scaling function and exponent. 
The latter is fixed by recalling that the Einstein relation 
connecting the growth of the characteristic length and 
the total transmitted intensity as a function of the linear 
size of the system [2^ |3l| gives I{L) ^ L^~^. On the 
other hand, plugging Eqs. (HI]) and into Eq. 1^ 

gives 

liL) = J T{t,L)dt = L'^ j f{L/l{t))dt - (23) 

which yields 77 = 1 — 2z. 

The transmission profile Tp {r, L) defined in Eq. ([8|) , 
does not depend on time, which was integrated over, but 
on two distances, r and L. The relevant scaling function 
is expected to depend only on the ratio of these distances, 



%{r,L)^L'^f{r/L) 



(24) 



Of course / and 77 are different from those appearing in 
Eqs. (|20p and (|22p . Since the integration of 7^(r, i) over 
all possible distances r gives once again the total trans- 
mitted intensity defined in Eq. (jlOp , the above-mentioned 
Einstein relation fixes the exponent to rj ~ —z. 



V. SIMULATIONS AND RESULTS 

In the present section we provide numerical evidence 
corroborating the scaling hypothesis and the single-long- 
jump picture illustrated in the previous section. An effec- 
tive encoding of the topology of the structure in a matrix 
of size g X L, where g and L are the generation of the 
carpet and the length of its side, respectively, allowed us 
to address L's up to a few million sites. In fact, after 
a straightforward generalization, the matrix relevant to 
a given carpet can be used to describe its d-dimensional 
counterpart with no additional memory demand. We will 
present some results for three dimensional structures as 
well. 

We start by focusing on the local probability density 
Pj{r,t). After checking that, as expected, the asymp- 
totic behaviour of this quantity does not depend on the 
starting site j, we collected extensive data on Po{r,t), 
where O denotes one of the vertices of the structure under 
examination. Figure |4] demonstrates that Pq {r, t) satis- 
fies the scaling hypothesis in Eq. (fTT|) . with a suitable 
choice of the exponent governing the asymptotic growth 
of the characteristic length £(t), as dictated by Eq. p^ . 
The expected steep decay of the scaling function f{x) is 
also apparent from Fig. 21 while the inset in the same 
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FIG. 4: (Color online) Scaling of the local probability for 
the (thin) SC with Ur = 8 and Uu ~ 3, corresponding to 
a ~ 0.8928 (head-on dynamics). The inset provides a check 
of Eq. (|13l) for some values of p. A fit of the data according 
to Eq. dni) gives z-^ « 0.553. 



figure demonstrates the correctness of our surmises in 
Eqs. (fT5|) and (fT^ . For ID systems, analytic arguments 
show that the exponent appearing in Eq. (|12p depends 
on the parameters n„ and determining the topology of 
the structure through the parameter a defined in Eq. ([3]), 
as z = 1 + a [1^ . Fig. [S] shows the dependence of z on 
a as obtained from fits of Eq. (|13p for two and three di- 
mensional structures. For the former we employed both 
dynamics, whereas for the latter only the head-on dynam- 
ics was considered. Notice that we examined structures 
characterized by different values of n„ and Ur, which com- 
bine to give very similar values of a, and that the cor- 
responding fitted values of z exhibit a very satisfactory 
agreement. This is rather convincing evidence that, as 
in the ID case, the exponent z depends on a rather than 
on Till and separately. Also, our data shows that the 
anomalous growth of £{t) - i.e. z < 2 - is present only 
for a < 1, whereas for a > 1 the growth of the charac- 
tcric length is consistent with z = 2. That is, a < 1 and 
a > 1 correspond to superdiffusive and standard diffusive 
behaviour, respectively. This applies independent of the 
system dimensionality and chosen dynamics. Finally, we 
notice that, for the head-on dynamics on 2D structures 
the dynamic exponent seems to be a very simple func- 
tion of the Levy parameter, z = 2/(2 — a), which might 
indicate that an analytic formulation is possible as in the 
ID case. 

We now turn to the properties of the average prob- 
ability in Eq. ([6|) . Figure [6] shows that the scaling re- 
lation in Eq. pT|) applies for sufficiently small values 
of r/i{t), where £{t) is the characteristic length dictat- 
ing the asymptotic behaviour of the local probabilities 



FIG. 5: (Color online) Dynamic exponent governing the 
asymptotic behaviour of the characteristic length in Eq. (|12p 
for several quasicrystals and jump schemes. Head-on dy- 
namics on 2D structures (circles); fan-out dynamics on 2D 
structures (triangles); head-on dynamics on 3D structures 
(squares). The solid line is the analytically known result for 
the ID case, z — a -\- 1, [2^, shown for comparison. The 
dashed line is an ansatz for the head-on dynamics on 2D sys- 
tems, z = z/(2 — a). The dash-dot line is the result for 
"ballistic" motion. Note that we examined structures char- 
acterized by very different geometries yet corresponding to 
very similar values of the parameter a in Eq. Q. Specifi- 
cally we considered a — log 8/ log 3 — 1 f» log 40/ log 7 — 1 
in the two-dimensional case and a = log 3904/ log 20 — 2 « 
log 218/ log 7 — 2 in the three dimensional case. 



Pj{r,t). For large values of r/£{t) it is possible to recog- 
nize structures ultimately arising from the complex su- 
perposition of delta functions appearing in Eq. (|14|) . The 
long tails characterizing the distributions P(r, t) give rise 
to the anomalous behaviour in the average moments de- 
scribed by Eqs p8| - p9| . These results of the single-long- 
jump hypothesis are indeed confirmed by Fig. [71 showing 
the fitted values of the exponent C dictating the asymp- 
totic dependence of the moments on time, (r^) ^ t'''-P\ 
The linear dependence of C{p) on p and the change in 
slope predicted by Eqs (fTS jl - lfTQ)) is clearly recognizable. 

We conclude the present section by discussing fur- 
ther results confirming the general validity of scaling ap- 
proach. Figs. [8] and [9] refer to the escape-time probability 
for the head-on and fan-out dynamics, respectively. It is 
apparent that S{L,t) obeys Eq. (PO)) . where -q = —z and 
z depends on a according to the behavior described in 
Fig. El 

The same applies to transmission-related quantities, as 
it is clear from Figs. [TOl and [TTl which prove the validity 
of scaling Equations ([22]) and ([24]), respectively. We note 
that Fig. [TTl shows that at small r/L values the conver- 
gence of the transmission profiles is attained only for suf- 
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FIG. 6: (Color online) Scaling of the average probability for 
the (thin) SC with rir = 8 and Uu ~ 3, corresponding to 
a ~ 0.8928 {head-on dynamics). The characteristic length 
£{t) is exactly the same as that used in Fig. (H £{t) ^ 



ficiently thick samples. At first sight this could be related 
to the cusp- like shape observed in the experiments 
However, when plotted in logarithmic scale, the profiles 
turn out to be exponential for large r/L and bell-shaped 
at the center. This is the same shape characterizing the 
transmission profile for a homogeneous sample (see inset 
of Fig. Ilip . The discrepancy between the experimental 
profiles and our results is very likely to be ascribed to 
the fact that the experimental samples are much more 
inhomogeneous than our Levy quasicrystals. 
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FIG. 9: (Color online) Scaling of escape time probability. 
Left: Q = log 24/ log 7 - 1. Right: a = log 16/ log 5 - 1. 
The numeric data have been obtained by adopting the fan- 
out dynamics described in Fig. [3l 



1 2 3 
P 



2 4 6 8 

P 




10" 



10" 



10" 




10" 



FIG. 7: (Color online) Exponent governing the asymptotic 
behaviour {r'') = f^'*"' of the moments of P{r,t), Eq. 0. 
Left: SC with a = log 8/ log 3 - 1; dotted line: ( = z^^p; 
dashed line: C — P ~ z~^{z — l)a; symbols: data fits. Right: 
FC with a = log 16/ log 3 — 1; dotted line: = p/2; dashed 
line: ^ = p — a -\- 1/2; symbols: data fits. 



FIG. 10: (Color online) Scaling of time-resolved transmis- 
sion probability. Left: a = log 24/ log 7— 1. Right: a — 
log 16/ log 5 — 1. The numeric data have been obtained by 
adopting the fan-out dynamics described in Fig. [3] 
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FIG. 11: (Color online) Scaling of the transmission profiles for 
the standard SC, a = log8/log3 — 1, according to Eq. (|24p . 
The inset shows a comparison between the profile for L = S*" 
(black thin line) and that for a homogeneous sample (gray 
thick line). The numeric data have been obtained by adopting 
the second dynamics described in Fig. [3] 



VI. CONCLUSIONS 

In this paper wc analyze Levy walks whose step-length 
distribution is characterized by correlations arising from 
the structural properties of the Levy quasicrystal through 
which the process takes place. One of the main outcomes 
of our work is that this correlation has a non-trivial in- 
fluence over all of the properties of Levy walks we have 
addressed. The difference between correlated and un- 
correlated Levy walks are not a peculiarity of the one- 
dimensional case examined in Refs. [26l - [28| . but persists 
in two- and three dimensions, at least as long as the corre- 
lations arise from the self-similar inhomogeneity pattern 
of the Levy quasicrystals hereby considered. The scaling 
analysis we develop shows that the dynamic exponent z 
governing the asymptotic behaviour of the characteris- 
tic length of the local probability Pj{r,t) has a wider 
scope than the mean-square displacement usually em- 
ployed in the description the process of (super-)difFusion 
through the sample. Our single-long-jump ansatz pro- 
vides a very satisfactory quantitative description of the 
strongly anomalous asymptotic behaviour of the process 
emerging after averaging over all possible starting points. 
In order to make a more direct contact with experiments 



we extend our scaling analyis to escape-time probabili- 
ties, time-resolved transmission probabilities and trans- 
mission profiles. Our scaling picture is indeed confirmed 
by extensive numeric simulations covering a significant 
range of sample sizes and heavy-tailed jump distribu- 
tions. 

A further important message from our results is that 
Levy quasicrystals are an ideal testbed for reproducing 
the superdiffusive features induced by structural proper- 
ties of the sample. In this respect the question naturally 
arises as to what extent the above scenario is robust to 
the introduction of disorder typical of the highly inho- 
mogeneous structures usually invoked for the realization 
of correlated Levy walks |10l-[l4l|. Specifically, the com- 
parison of our transmission profiles and those reported 
in Ref. [3] suggests that Levy quasicrystals are much 
more homogeneous than the Levy glass realized at Eu- 
ropean Laboratory for Non-Linear Spectroscopy. We are 
currently working at exending our analysis to structures 
characterized by disorder and/or a higher degree of in- 
homogeneity than Levy quasicrystals [331 ■ 

While we deem these issues of utter interest, we think 
that Levy quasicrystals should not be regarded as a sim- 
ple toy model, but as an alternate, highly controllable 
route to observe correlated Levy walks in engineered ma- 
terials. In this perspective, the "cleanliness" of the sam- 
ples suggested by the bell-like shape of the transmission 
profile in Fig. [TT] could even prove to be a desirable fea- 
ture. Indeed, in the absence of the strong inhomogeneity 
characterizing most of the relevant physical systems, the 
superdiffusive character of Levy walks could emerge more 
cleanly from experimental measurements. For instance 
a data-collapse of measured time-resolved transmission 
probabilities like in Fig. [TU] would provide a measurement 
of the dynamical exponent characterizing the asymptotic 
behaviour of the process. 

The comparison of measurements on Levy glass and 
Levy quasicrystals would then highlight the role of strong 
inhomogeneity in quenched Levy walks. 
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